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Scalar order: possible candidate for order parameters in skutterudites 

Annamaria Kiss* and Yoshio Kuramoto ' 
Department of Physics, Tohoku University, Sendai, 980-8578 

Phenomenological Landau analysis shows that the properties of ordered phases in some skutterudites 
are consistently accounted for by a scalar order parameter which preserves the cubic symmetry, even in 
the ordered phase. A universal value is found for the anisotropy ratio of the transition temperature in 
a magnetic field, homogeneous magnetization, and induced staggered magnetization. The difference in 
magnetic behavior between PrFe4Pi2 and PrRu4Pi2 near their phase transitions is explained within a single 
framework. For the low-field phase of PrFe 4 P| 2 , the scalar order with the T tg symmetry can explain (i) the 
absence of field induced dipoles perpendicular to the magnetic field, (ii) isotropic magnetic susceptibility 
in the ordered phase, (iii) the field angle dependence of the transition temperature, and (iv) the splitting 
pattern of the 31 P nuclear magnetic resonance (NMR) spectra. It is proposed how the order parameter in 
SmRu 4 P|2 is identified by NMR analysis of a single crystal. 
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Rare-earth filled skutterudites have been attracting consid- 
erable attention from both experimental and theoretical sides 
because of their intriguing behaviors. Among them, Pr-based 
compound PrFe 4 Pi2 shows a phase transition at Tq = 6.5K, 
which can be seen as sharp anomaly in the magnetic suscepti- 
bility. 1 In the ordered phase staggered dipoles are induced by 
magnetic field, 2 which suggests that the order parameter does 
not break the time-reversal symmetry. PrRu 4 Pi2 has a metal- 
insulator phase transition at Tmi = 65K, and its crystalline 
electric field (CEF) states show drastic change below 7mi, 3 
which seems to be described by antiferro-type order of hex- 
adecapole moments with Fi f symmetry. 4 S111RU4P12 has also 
a metal-insulator phase transition at 7mi = 16. 5K. The na- 
ture of the order parameter in phase II is not clear until know. 
An octupolar order with the Fs u symmetry, which breaks the 
time-reversal symmetry, has been proposed for this phase. 5 In 
fact, nonzero internal field was observed below Tmi by recent 
yuSR experiment. 6 

In this paper we propose that both PrFe4Pi2 and PrRu 4 Pi2 
have a scalar-type order parameter with the Ti f symmetry. 
We show by phenomenological analysis that the scalar order 
model explains the main properties of PrFe4Pi2 including the 
NMR results: (i) the absence of field induced dipoles perpen- 
dicular to the magnetic field, (ii) isotropic magnetic suscepti- 
bility in the ordered phase, (iii) the field angle dependence of 
the transition temperature, and (iv) the splitting pattern of the 
31 P NMR spectra. As a first step, we concentrate on behaviors 
at low magnetic fields. For SmRu 4 Pi2, in view of its nearly 
isotropic behavior in the ordered phase, we propose another 
candidate of the octupole order of the r nz type, which trans- 
forms as a pseudo-scalar with the Ti,, symmetry. It is proposed 
how the order parameter in SmRu 4 Pi2 is identified by NMR 
for the single crystal. 

Up to the present, the order parameter in PrFe4Pi2 has 
widely been considered as an antiferro-quadrupolar(AFQ) or- 
der of T3 moments. However, this AFQ model fails to account 
for the isotropic susceptibility in the ordered phase, for exam- 
ple. Furthermore, with static F3 quadrupoles it is difficult to 
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explain why the field induced staggered dipoles are always 
parallel to the field direction, as indicated by neutron diffrac- 
tion 2 and NMR. 7 ' 8 

A scalar-type order can be of two different kinds: one (Ti u ) 
breaks, while the other (F 1? ) does not break the time-reversal 
symmetry. On the other hand, both of them preserve the cubic 
symmetry even in the ordered phase. Therefore, the Landau- 
type expansion of the free energy T contains cubic invariants 
composed by the magnetic field components. Around a sec- 
ond order phase transition we expand f up to fourth order in 
the order parameter and magnetic field as 



1 



1 



rOfci, H) = T Q (H) + -a s [T - T c (H)]ifr 2 Q + -b s ifr 4 Q , (1) 

where i^q is the staggered component of the scalar order 
parameter with ordering vector Q = [1,0,0], and H = 
H(h x , h y , h z ) is the external magnetic field. The transition tem- 
perature in magnetic field has a dependence: 

T C (H) = T + ^t 2 H 2 + i(r 4 + t 4a h 4 )H 4 , (2) 

where £2, h, f 4fl are expansion coefficients, and /z 4 =h\ + h 4 + 
h 4 -3/5. This invariant /z 4 is common in the case of cubic (O/,) 
and the tetrahedral (7 1 /,) point group symmetries. Therefore, 
our present treatment is valid for both cases. The first part 
Tq{H) has a field dependence similar to eq.(2). 

The anisotropy in eq.(2) is independent of the microscopic 
details of the scalar order. Using /i 4 [100] = 2/5,/z 4 [110] = 
-1/10,/! 4 [111] = -4/15 we obtain the ratio 

r e [Ooi]-r c [iii] 



= 4. 



(3) 



r e [iio]-r c [iii] 

This relation should hold as long as the magnetic field 
is weak enough when the scalar order emerges. Figure 1 
shows the transition temperature in PrFe 4 Pi2 measured as a 
function of field angle 9 which is defined by {h x ,h y ,h z ) - 
(cos sin 0, sin</>sin#, cos#). The anisotropy of eq.(2) with 
f 4fl < provides excellent fit to the observed T c . This result 
strongly suggests the scalar order in this compound. 

Equation (1) also shows that the magnetization along three 
principal axes should have the anisotropy ratio: 

(Afioo-M m )/(Afii -M m ) = 4, (4) 
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or, equivalently, (M m - Afioo)/(Afno - Afioo) = 4/3. The ra- 
tio given by eq.(4) holds both in the paramagnetic and ordered 
phase as long as the magnetic field is weak enough. Because 
of the coupling with ^q, however, the weight of anisotropic 
part should change at the transition temperature. It is possi- 
ble that anisotropy is reversed by keeping the ratio given by 
eq.(4). Experimentally, the reversed anisotropy is indeed ob- 
served at T = 0.3K. 9 

We now discuss anisotropy of the staggered magnetization 
niQ with the order parameter F lg . The change of the free en- 
ergy due to niQ is given by 

T(<{fQ,mQ,H)-T(t{fQ,H) 

= ]ja m m\ + <AQ»i_Q \c\H + (c 3 + C3 a h 4 )H 3 j , (5) 

where T(i//q,H) is given by eq.(l). In the presence of 
with symmetry Fi g , itiq is induced by external magnetic field. 
Terms with coefficient c, in expression (5) come from the in- 
variant Fig(Q) ® r4„(-Q) ® r4„(0), which requires that the 
induced dipoles are always parallel or anti-parallel to the field 
direction. From the condition df/dmQ = we obtain 

m Q = -— [ciH + (c 3 + c 3a h 4 )H 3 j if/ Q , (6) 

which shows explicitly that niQ develops in magnetic field 
when the order parameter i//q is non-zero. In the linear order 
of magnetic field, there is no anisotropy in the induced stag- 
gered dipoles. In the third order, the anisotropy ratio is again 
given by 

(m Q [100] - m Q [l 1 1])/(otq[1 10] - m Q [l 1 1]) = 4. (7) 

At present, there is no experimental information for the above 
ratio, since neutron scattering has been done only for mQ[100] 
and »iq[1 10]. 2 The same anistropy ratio given by eqs.(3), (4) 
and (7) comes from the property of h 4 alone, and is a universal 
feature for the scalar order. 

Let us explain how the difference between PrFe4Pi2 and 
PrRu4Pi2 can be interpreted in a phenomenological frame- 
work. The most obvious difference appears in the magnetic 
susceptibility; PrFe4Pi2 shows a sharp peak at the transition, 1 



while no conspicuous anomaly is seen in PrRu4Pi2- In consid- 
ering the magnetic susceptibility around the zero-field phase 
transition temperature To, it is more convenient to use the 
Gibbs potential which is obtained from the free energy T 
by Legendre transformation. Namely we obtain using the ho- 
mogenous dipole moment M = mo, 

= T^ Q ,H = 0) + l -a f {T - T F )M 2 + ^Aifr 2 Q M 2 , (8) 

where Tity^H = 0) is given by eq.(l) with H = 0, and the 
last term represents the coupling between the scalar order and 
the magnetization. The term with a f becomes important if the 
ferromagnetic instability is close to the scalar phase transi- 
tion: Tp < Tq, as in the case of PrFe4Pi2. Indeed, a positive 
Curie-Weiss temperature Tp « 3.5K has been found in exper- 
iments. 1 To the contrary, PrRu4Pi2 seems far from the ferro- 
magnetic instability. In this case, M does not have a significant 
amplitude, and the coupling term with A is less significant. 

Around the transition temperature To, the magnetic suscep- 
tibility x+ for T > To and^_ for T < To can be expressed as 

X - X =a f (T -T F ), (9) 
XZ 1 = (a f - Aa s /b s )T - a f T F + Aa s T Q /b,, (10) 

where both x± are isotropic and follow the Curie- Weiss law. 
If the coupling with magnetic moment is strong enough, we 
obtain Aa s /b s > a f , which leads to a peak in the susceptibility 
at Tq. Figure 2 shows comparison between theory and experi- 
ment for PrFe4Pi2- The agreement is excellent with the choice 
of Aa s /(afb s ) = 2.8. In consistency with the sharp peak in the 
susceptibility, the phase boundary is suppressed appreciably 
by magnetic field in PrFe4Pi2. 

In the opposite limit of negligible A, xZ l given by eq.(10) 
reduces to x~+ ■ Hence there is no change at the scalar 
phase transition. This limit seems to explain the situation in 
PrRu 4 Pi 2 . Accordingly, the transition temperature is hardly 
affected by magnetic field. 10 

We proceed to analyze 31 P NMR spectra in skutterudites 
in terms of the scalar order. We discuss mainly the case of 
PrFe4Pi2, but also touch on SmRu4Pi2 where the pseudo- 
scalar Fiu is a candidate of the order parameter. Around each 
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Pr ion at position (0,0,0), there are six P positions r i( 2) = 
(0, v, ±u), T3 ( 4) = (+u, 0, v) and r 5( 6) = (v, ±u, 0), which 
are crystallographically equivalent. With finite magnetic field, 
these positions are no longer equivalent, and splitting of NMR 
lines takes place. The splitting is of purely magnetic origin, 
since 31 P ions have no quadrupolar moment with / = 1/2. 
Therefore, interactions are only between the P nuclear spin I 
and the dipole J and octupole T^, T xyz moments of a Pr ion. 

Taking a representative pair at r 3 and r 4 , we write down 
the hyperfine interaction characterized by energies for the 

dipoles and e^J for octupoles with k, I distinguishing indepen- 
dent components. 11 The invariant form of the interaction con- 
tains terms such as I X J Z and I y T xyz since there is only a mirror 
symmetry against the xz-plane for the pair. The hyperfine in- 
teraction is given by 

ff M (3,4) = I x [ef\j x ± e%J z + e£>7* ± e^Tf] 



+i z [ep z ± ep x + ^Jrf ± 4°jr?] , (11) 

where the negative sign corresponds to P ion at position r 4 . 
The interaction for other pairs can be obtained from eq.(ll) 
by using proper rotational operations. 

Let us discuss first the case of PrFe4Pi2. In Tf, symmetry, 
the dipole moment J corresponds to rJJ and the octupole mo- 
ment T* 3 to T®, but they are mixed due to the lower symmetry. 
Therefore, both are induced by the external magnetic field. In 
the disordered phase (<//q = 0), the homogenous moments are 
induced as m (= Jo) = MH/|H| and = TH/|H|, where M 
and T are magnitudes at a given temperature and magnetic 
field. They cause the splitting of the P NMR line depend- 
ing on the field direction. Since the external magnetic field is 
much larger than the hyperfine field, we assume (I x , /,,, I z ) oc 
(h x , h y , h z ). In the case of UFe4Pi2, classical approximation of 
the dipolar field shows good qualitative agreement with the 
measured NMR spectra. 12 However in PrFe4Pi2, we checked 
that the approximation leads to a large deviation from the 
measured result. 7 Hence we fix the parameters in eq.(ll) in 
a phenomenological manner. We define the hyperfine field hhf 
so that //hf = y n Ihhf = I(f - /o), where y n is the nuclear 
gyromagnetic ratio of 31 P, / is the resonance frequency and 
/o is the zero shift. We write /ihf(l,2) = g\, /ihf(3,4) = g 2 
and Ahf(5,6) = g 3 for H || (001), and Ml, 3,5) = k x for 
H || (111). We obtain 

Ml, 2) = g 3 h 2 x + g 2 h 2 y + gilt] + g 4 h y h z , 

M3, 4) = g\h\ + g 3 h 2 y + g 2 h 2 z + g 4 h z h x , 

M5, 6) = g 2 h 2 x + gih 2 + g 3 h 2 + g 4 h x h y , (12) 

where g 4 = 3k\ - g\ - g 2 - g 3 and parameters gj (I < i < 4) 
are linear combination of ef j (a = d, o) times M or T. We 
determine the four parameters g, by fitting to the observed 
three lines for H || (001) and the three degenerate ones from 
PI, P3, P5 for H || (111). Figure 3 shows our results for the 
fitting in the disordered phase. The values for parameters gi 
are summarized in Table I. The result for H || (110) is a con- 
sequence of the form of the hyperfine interaction given by 
eq.(ll), which is independent of the microscopic details of 
the model. The spectrum computed for H || (1 10) is found to 
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Fig. 3. Theoretical fitting (six lines) to experimental 31 P NMR results 7 
(crosses) in PrFe4Pi2 at magnetic field H = 4T and temperature T = 50K 
as a function of field direction with 4> = n/4. The numbers indicate the P 
positions. 



Table I. Choice for the parameters to describe the measured NMR results. 



gl 


g2 


g3 


h 


5.2 mT 


10.7 mT 


0.2 mT 


6.3 mT 


a\ = a-i 


a 2 


c\ 




1.2 x icr 2 


0.8 x 10~ 2 


i.i4 x ur 1 


o.99 x ur 1 


dt 


di 


d 3 


d 4 


1.0 x 10^ 


1.2 x 10~ 2 


1.1 x 10^ 


o.9 x ur 2 



be in reasonable agreement with experimental results. Since 
eq.(ll) includes only nearest neighbor interaction between 
the Pr and P ions, the slight deviation between theory and 
experiment is ascribed to effects of distant Pr-P pairs. 

Let us now consider the case of i^q + with the symme- 
try ri g . We restrict our discussion to the case of low magnetic 
fields, and expand the quantities in linear order of magnetic 
field. Thus we write the staggered dipoles and octupoles in- 
duced by the magnetic field as diq = K\ ,</</>qH and Tq = 
K ho \j,QK. As a result, Pr positions (0, 0, 0) and (1/2, 1 /2, 1 /2) 
become inequivalent, and extra splitting of the NMR lines de- 
velops below T c . In magnetic field along (001), for example, 
the extra splittings of the three main lines are described by 

AMI, 2) = (Aq (Kite™ + K Uo ef})H = a x H, 
AM3.4) = <AQ {K l4 ef x + K Uo ef})H = a 2 H, 

AM5, 6) = iAq (K i4 e™ + K ho e^ H = a 3 H. 

We determine magnitudes of parameters a, (i = 1 , 2, 3) so as 
to reproduce the corresponding experimental results. 13 Figure 
4 shows the result of fitting together with experimental results. 
It is obvious that the experimental results deviate from the lin- 
ear behavior for H > IT The different field dependence for 
the splittings AMT, 2)(« Mm(5, 6)) and A/i hf (3, 4) is intrigu- 
ing. It seems hard to understand the difference without con- 
sideration of induced staggered octupoles. The analysis with 
non-linear effects of magnetic field is rather complicated, and 
will be presented in a separate publication. 

For H || (1 1 1), we define the parameters c\ for A/ihf(l, 3, 5) 
and c 2 for A/zm(2,4, 6) in a way analogous to a,-. The sim- 
ple splitting of lines (1,3,5) and (2, 4, 6) in the ordered phase 
can be explained with scalar order, but not with quadrupo- 
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Fig. 4. Magnetic field dependence in linear order in H (solid line) of the 
extra line splittings for field direction (001) in PrFe4Pi2. Symbols are the 
measured result at T = 2K 8 and numbers indicates the P positions. The 
parameters used are shown in Table I. 

lar order. The splittings A/ihf(l,3,5) and A/ihf (2, 4, 6) have 
almost identical field dependence with a tiny deviation at 
low fields. 7 We fix the parameter c\ to reproduce the ex- 
perimental result, and then parameter is determined as 
C2 = (2ai+2a2 + 2fl3-3ci)/3.ForH || (110), we define d\ for 
Mhf(l,2), di for A/i hf (3,4), J 3 for A/i hf (5), andc/ 4 for A/z hf (6). 
It turns out that all di can be fixed by the parameters aj and c\ 
as d\ = (ci2 + aj)/2, di = (a\ + ai)j2, dj, = (3ci - «3)/2 and 
di = (2a i + 2«2 + «3 - 3ci)/2. Namely, experimental results 
along (110) should be reproduced without further adjustable 
parameters provided eq.(l 1) applies to the actual system. Ex- 
perimental values with superscript e are c e 2 = 1 .0 x 10~ 2 , d\ ~ 
d\ = 0.85xl0~ 2 , d\ = 1.2 x 10~ 2 , d\ = 1.1 x 10~ 2 [7], which 
show reasonable agreement with corresponding theoretical 
values shown in Table I. We expect qualitatively similar pat- 
tern for the NMR spectra in the ordered phase of PrRu4Pi2, 
since its order parameter should also be a scalar. 

We now discuss briefly the expected splitting pattern of the 
NMR spectra in the case of T xyz octupolar order. When the 
magnetic field is applied along (001), the hyperfine interaction 
for different pairs is given as 

H hs (\,2) = e^ l I z J z , 
Hhf(3,4) = 4*/ Z 7 z , 

H M (5, 6) = e ( 2 %J z ± e ( 2 %T xyz . (13) 

The T xyz octupolar order causes the splitting of lines (5,6). 
With the octupole order parameter we obtain A/ihf(5) = 
A/ihf(6) oc e 2 ° 2 iJ/Q. The splitting occurs also for the octupole 
ordering vector q = because &m(5) - hhi (6) oc e£ 2 >^o + in 
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this case. Then four NMR lines appear in the ordered phases 
with q = Q as well as q = 0. However, in the case of field 
along (111), doubling of each line occurs with q = Q, while 
no extra splitting is expected with ordering vector q = 0. The 
former splitting pattern is similar to the case of PrFe4Pi2. 

In this paper we have considered the characteristic features 
of scalar orders with symmetries Y\ g and T\ u in weak mag- 
netic fields. We have found the universal anisotropy ratio in 
weak magnetic field by phenomenological Landau-type anal- 
ysis. We conclude that the scalar order scenario with sym- 
metry Y\ g can explain the known properties of PrFe4Pi2 con- 
sistently. We have shown that the splitting pattern of the 31 P 
NMR spectra in the disordered and ordered phases can also 
be described within this framework. We have predicted the 
31 P NMR spectra in the case of T xyz octupole order for fields 
along (001) and (111). This octupole moment transforms as a 
pseudo-scalar with the T\ u symmetry and can be a good can- 
didate for the order parameter in SmRu 4 Pi 2 below the metal- 
insulator transition. 
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